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基于 线性 与 非 线性 特征 融合 的 J 波 自 动 识别 - 
组 凡 铭 ， 李 灯 熬 '"， 赵 菊 敏 


(太原 理工 大 学 信息 与 计算 机 学 院 ,， 山西 晋 中 030600) 


摘 要 : 临床 研究 表明 , J] 波 可 作为 一 些 心脏 疾病 的 高 危 预警 指标 。 针 对 当前 医生 诊断 J 波 仅 通过 经 验 进行 识别 ， 易 造 
成 误诊 的 问题 ， 从 信号 处 理 角度 ， 提 出 了 一 种 波 自动 识别 方法 。 该 方法 通过 提取 心 电 数 据 极点 对 称 模 态 分 解 后 的 能 
量 特 征 与 高 阶 累 计量 特征 ， 融 合 线性 与 非 线 性 特征 ， 并 采用 主 成 分 分 析 进 行 特征 降 维 ， 最 后 利用 人 工蜂 群 算法 优化 后 
的 支持 向 量 机 进行 分 类 ， 实 现 丁 波 的 自动 识别 。 对 比 实 验 结 果 ， 提 出 的 方法 平均 准确 率 为 97.3%， 可 有 效 地 识别 J 波 。 
关键 词 : J 波 ; 极点 对 称 模 态 分 解 ; 高 阶 累计 量 ; 人 工蜂 群 算 法 ; 支持 向 量 机 
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Automatic recognition of J wave based on combination of linear and nonlinear features 


Wu Fanming, Li Dengao', Zhao Jumin 


(College of Information & Computer, Taiyuan University of Technology, Jinzhong 030600, China) 


Abstract: Clinical studies have shown that J wave can be used as a high risk early warning index of some heart diseases. In 
view of the shortcomings of J wave diagnosis only by the clinicians’ experiences, which can easily lead to misdiagnosis. From 
the perspective of signal processing, this paper proposed a J wave automatic identification method. This method extracted the 


energy features of electrocardiogram data after the extreme-point symmetric mode decomposition and higher order statistics, 


combined linear and nonlinear features, and adopted the principal component analysis to reduce the dimension of the features. 


Finally, it used support vector machine optimized by artificial bee colony algorithm to realize automatic identification of J wave. 


The experimental results Show that the average accuracy rate of identifying J wave is 97.3%, which can effectively identify the 
J wave. 
Key words: J wave; extreme-point symmetric mode decomposition; higher-order statistics; artificial bee colony algorithm; 


support vector machine 


蜂 群 算法 (artificial bee colony，ABC) 改进 后 的 支持 向 量 机 
(support vector machine，SVM) 对 心 电 信 号 进行 分 类 ， 从 而 实现 
J 波 的 自动 识别 。 医生 依据 ] 波 自动 识别 的 结果 , 并 结合 患者 的 
临床 状况 ， 提 出 最 佳 诊 断 和 治疗 方案 。 


A 


丁点 是 心 电 信 号 上 QRS 波 与 ST 段 的 结合 点 ， 表 示 心 脏 队 
极 的 结束 、 复 极 的 开始 帆 。 若 丁点 从 基线 移 位 ， 呈 尖峰 状 、 驼 
状 或 顿挫 状 , 且 振 幅 >=0.1mv、 时 限 >=20ms, 该 偏 移 波 即 为 了 波 


针 


ZI 


1 ”数据 准备 及 预 处 理 


[2,3], 
相关 研究 表明 ， 当 J 波幅 度 超过 0.1mv 且 伴 有 ST 段 变 化 本 文 从 山西 大 医院 心 电 数据 库 中 获得 了 大 量 含 J 波 的 纸 质 
时 ， 极 易 引 发 致命 心肌 梗塞 、 恶 性 心律 失常 ， 甚 至 心脏 性 独 。” 心 电 数 据 。 使 用 高 分 辩 率 扫描 仪 将 纸 质心 电 数据 转换 为 数字 图 


i 
攻 


死 1。 当 前 临床 上 医生 通过 观察 心电图 上 了 波 的 大 小 、 时 限 并  ” 片 ， 并 采用 图 像 处 理 技术 ,对 转换 后 的 数字 图 片 进行 Hough 变 
结合 患者 的 临床 表现 进行 鉴别 诊断 ， 依 据 主观 经 验 来 进行 J] 波 换 、K-means 聚 类 、 曲 线 拟 合 等 , 实现 纸 质心 电 数 据 的 数字 化 ， 


识别 ， 容 易 造 成 误诊 和 漏诊 ， 因 此 实现 了 波 的 自动 识别 ， 可 及 ”具体 步骤 如 下 : 

时 对 本 波 进行 诊断 ， 有 效 降低 了 波 的 致死 率 。 a) 倾斜 校正 。 采 用 霍 夫 (Hough) 变换 检测 心电图 纸 中 背 
为 了 实现 了 波 的 自动 识别 ， 本 文 从 信号 处 理 角 度 对 J 波 研 。 景 网 格 直 线 ,通过 计算 斜率 得 到 倾斜 角度 ,进而 进行 旋转 校正 。 

究 ， 通 过 提取 心 电 数据 的 线性 特征 和 非 线性 特征 ， 并 通过 人 工 b) 像素 聚 类 。 同 为 心 电 数据 或 同 为 背景 网 格 像素 具有 相似 
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性 ， 本 位 采用 K-means 每 


案 类 算法 ， 分 离 黑 色 的 心 


的 背景 网 格 。 


c) 移 除 干扰 。 由 于 纸 质心 电 数据 存在 导 联 标志 、 
等 干扰 像素 点 ， 所 以 从 二 值 图 像 中 选择 联通 
进行 移 除 。 


有 数据 和 红色 


基准 电压 


d) 


线 拟 合 。 读 取 二 值 图 像 中 心 


标 ， 并 对 缺失 


的 像素 点 


化 后 的 完整 心 
以 250 Hz 的 采样 
常 心 电 数据 和 40 个 


数据。 


率 从 数字 化 后 的 心 


根据 原 记 录 分 别 标记 为 类 型 N 和 类 型 J。 


于 心 电 数 据 的 频谱 主要 


由 于 外 界 因素 干扰 ， 心 电 数据 中 通常 包含 多 种 噪声 四。 由 


数据 的 各 像素 点 位 置 华 
进行 插值 。 通 过 拟 合 像素 点 ， 获 得 数字 


有 数据 中 选取 40 个 正 
含 J 波 心 电 数 据 ， 每 个 记录 持续 2 min， 并 


面积 较 小 的 像素 点 


l 


扰 。 


2 ”J 波 自动 识别 


本 文通 


过 对 预 处 理 后 的 心 电 数 据 提取 线 ; 


性 特征 


中 在 1 Hz 一 40 Hz 间 ， 本 文通 过 截止 
Hz 和 40 Hz 的 二 阶 巴 特 沃 斯 滤波 器 去 除 心 电 数 据 中 
的 基线 漂移 和 工 频 干 


与 非 线性 特 


疹 YZ 


全 ， 沐 上 


decomposition, ESMD ) 提取 线性 


特征 高 阶 


行 降 


特征 , 融合 心 


极点 对 称 模 态 分 解 (extreme-point symmetric mode 
有 数据 的 非 线 性 


累计 量 (higher-order statistics，HOS ) 并 对 融合 特征 进 


住处 理 , 将 采用 主 成 分 分 析 (principal component analysis， 


PCA ) 降 旨 


别 。 该 算 ; 


法 的 整体 框图 如 图 


1 所 示 。 


人 后 的 特征 输入 到 ABC 算法 优化 后 的 SVM 中 进行 识 


线性 特征 
ESMD 


PCA 降 维 一 一 站 ABC-SVM 


类 别 


2.1 


的 本 征 模 


数据 准备 
及 预 处 理 


线性 特征 
本 文采 / 


非 线性 特征 
HOS 


类 别 N 


pa 
站 
= 
[FHl 

对 


ESMD 来 提取 心 电 数据 的 线性 特征 ， 
态 函 数 (intrinsic mode function, IMF ) 的 能 量 特 行 


目前 对 心 


empirical mode decomposition,EMD ) 两 种 方法 。 


变换 在 心 


自 适 应 的 ， 


其 性 能 的 好 坏 严重 依赖 于 小 波 函 数 的 的 选择 ， 


初始 选 定 
达到 最 佳 
据 信 号 自 


小 波 函 数 ， 全 局 信号 必须 都 使 用 相同 的 小 波 基 ， 


提取 分 解 后 


征 [8]。 


数据 进行 时 频 分 析 的 方法 主要 有 小 波 变 换 和 经 验 模 


小 波 


数据 识别 分 类 有 很 好 的 分 析 效 果 ， 但 小 波 变 换 是 非 


一 旦 
很 难 


分 析 效 果 ; EMD 是 一 种 自 适应 的 信号 处 理 方法 ， 


ChinaXiv 盒 作 期 和 


毋 凡 铭 ， 二 线性 YR, 了 丁 ; 


ESMD 具体 分 解 过 程 如 下 : 

a) 找 出 原始 心 电 数据 x(?) 中 的 所 有 极 值 点 , 包括 极 大 值 点 
与 极 小 值 点 ， 记 做 PU <i<n)。 

b) 用 线段 连接 所 有 临近 的 ,并 将 连接 的 线段 的 中 点 标记 
为 TQ<izn-D ， 同 时 在 左右 两 端 添 补 边界 中 点 五 和 了 。 


c) 通 过 n+1 个 中 点 来 构建 s 条 内 插曲 线 世 ,…, 世 (s>D , 计 
算 其 平均 值 了 =(L+…+L)/s。 
d) 对 x(D) -工序 列 重复 上 述 三 个 步骤 ， 直 至 |L|<e (是 


允许 误差 ) 或 者 筛选 次 数 达 到 预 设 的 最 大 值 K ， 得 到 第 一 个 
IMF 分 量 c 。 

e) 对 剩余 序列 xD-c 重复 上 
r 为 单一 信号 或 不 再 大 于 预先 
IMF 分 量 6,G3,.….。 

f) 在 限定 区 间 [Kwin, Kwox] 内 改变 KK 值 ,重复 上 述 五 个 步 又 。 
然后 计算 序列 xD -7 的 方差 o? ,并 对 o/o。 和 天 进行 绘图 ( ov 
是 x(?) 的 标准 差 ) ， 在 图 中 找 出 o/o6 最 小 值 对 应 的 Ko ， 以 R 
作为 限制 条 件 再 次 重复 上 述 五 个 步骤 ， 最 后 分 解 至 只 剩 一 定数 
量 极点 的 剩余 项 r+ ， 就 是 原始 信号 x(〗) 的 自 适 应 全 局 均线 ， 亦 
即 趋 势 项 。 


四 个 步骤 ， 直 到 剩余 序列 
a 便 可 分 别 得 到 


经 过 分 解 ， 原 始 信 号 x(?) 可 表示 xD = > c+r ， 即 利用 


i=1 


系列 IMF 和 一 个 趋势 项 ~ 。 


ESMD 将 时 间 序 列 xD 分 解 成 一 


心 电 数据 对 应 的 方差 比率 随 最 大 筛选 次 数 的 变化 


104c 9 


方差 比率 (%) 


1 1 1 1 1 1 1 
0 5 10 15 20 25 30 35 40 


第 选 次 数 
图 2 心 电 数据 对 应 的 方差 比率 随 最 大 筛选 次 数 的 变化 


ESMD 分 解 的 一 大 改进 是 最 大 筛选 次 数 由 自 适 应 全 局 均线 


决定 。 分 解 前 需 设 定 剩余 极点 个 数 MM , ， 这 决定 了 趋势 项 7 的 


次 数 ， 其 形状 能 否 反映 数据 变化 趋势 全 赖 此 设 定 。 为 


其 根 


身 的 波动 规律 将 信号 分 解 为 从 低频 到 高 频 不 同时 间 尺 


度 分 量 ， 


ESMD 在 EMD 


j 内 部 极 
分 解 的 余 量 
为 整个 信 
次 数 ， 有 


Da 


于 心 电 信 号 的 研 
的 基础 上 ， 


更 加 适 | 究 分 析 。 


点 对 称 插值 法 蔡 代 EMD 中 的 三 阶 样 条 插值 ， 


里 三 赋 予 J 意义 ， 
号 的 “ 自 适 应 全 局 均线 ”， 
效 地 解决 了 EMD 中 的 模 态 混 堵 


问题 。 


通过 改变 外 包 络 线 的 插值 法 ， 


目 给 


通过 最 小 二 乘 思想 优化 余 量 ， 使 其 称 
从 而 确定 分 解 的 最 佳 


筛选 


满足 边界 性 插值 要 求 ， 需 使 MI, 之 4 。M , 如 果 选 择 过 大 ， 会 


降低 模 态 分 辨 率 ， 本 文 M , 取 4 个 。 通 过 对 G/a 入 进行 


绘图 为 图 


2。 在 图 中 找 出 CV an 最 小 值 对 应 的 玉 =36， 则 心 电 


信号 对 应 36 次 筛选 的 最 佳 分 解 结果 在 图 3 中 显示 , 共 分 解 为 7 
个 IMF 和 1 个 趋势 项 7 。 
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录用 稿 毋 凡 铭 ， 等 : 基于 线性 与 非 线性 人 oY /= 波 自 动 识 
本 心 电 数据 在 36 次 筛选 - 5 的 ESMD 分 解 结 寺 果 t 2 
-OU 瓦 = 之 cloj =12… 
| a 2 5 1 10= 

时 间 (s) 时 间 (s) 

i 0.5 0.5 g 2 
和 0 WV—h— VN 0) UAVASA VA 人 A 
三 三 民生 rl(t (2) 

多 1 2 3 a 1 2 3 . 和 2 (0 
时 间 (s) 时 间 (s) 

入 0.2 te 0.2 
EoMAAANV 信人 人 八 人 ~ c) 得 到 心 电 数据 的 能 量 向 量 : 
-0 -0.2 

0 1 2 3 0 于 2 3 
时 间 (s) 时 间 (s) 

0 ,~ 0.05 V= EB ,EFE, 已 |] (85 
om oN 仆仆 
Wi 1 2 > 0 1 2 3 

ye 对 向 量 VV 做 归 一 化 处 理 ， 得 到 
四 3 心 电 数 据 在 36 次 筛选 下 的 ESMD 分 解 结果 V'=[p,, p,,…, p,, p, | (4) 
根据 ECG 的 波形 特点 ， 观 察 图 3 可 以 看 出 每 一 个 IMF 所 代 中 
表 的 意义 : 前 两 个 分 量 c| ，C, 主要 代表 了 频率 最 高 的 QRS 波 
E=O,E +E), 

的 成 分 ; 从 C3 开始 有 P 波 开 始 分 解 ; 从 cj 开始 加 入 T 波 成 分 ; 四 

Cs 代表 P，QRS，T 波 低频 成 分 分 量 琶 加 ，Cc6, cy 代表 更 大 时 间 Pe 0 oe 

尺度 上 的 心脏 生理 调节 节律 。 从 图 4 中 可 以 看 出 趋势 项 + ， 即 

自 适应 全 局 均线 ， 很 好 地 反映 了 ECG 信号 的 变化 趋势 。 p, = 已 /已 ,及 , 表 示 IMF 和 7 的 能 量 在 心 电 数据 总 能 量 

心 电 数据 的 自 适 应 全 局 均线 中 的 比重 。 
一 一 心 电 数据 
~ ~ 趋势 项 r 
“| 能 量 向 量 V' 的 分 量 从 P 到 p, 代表 的 频率 逐渐 降低 ， 本 

Bl 文 共 提取 以 上 8 个 线性 特征 。 心 电 数据 类 型 N 和 类 型 ] 的 能 量 
和 向 量 分 布 对 比如 图 5 所 示 。 从 图 中 可 以 看 出 心 电 数据 类 型 N 能 
一 量 向 量 从 高 频 到 低频 变化 趋势 呈 降 低 趋 势 ， 能 量 主要 集中 在 高 


频 部 分 ， 心 电 数据 类 型 ] 能 量 向 量 能 量 最 高 出 现在 ps 。 


Mo 0 1 区 1 2 2 3 请 类 型 N 心 电 数 据 与 类 型 J 心 电 数 据 ESMD 能 量 向 量 图 对 比 
寺 间 (s) “人 
图 4 心 电 数据 的 自 适应 全 局 均线 3 
心 电 数据 通过 ESMD 分 解 后 得 到 一 系列 IMF 和 一 个 趋势 项 wd 只 
>、 本 
7 ， 计 算 IMF 和 的 能 量 ， 可 以 计算 得 到 能 量 向 量 。 计 算 过 程 于 。。| ek 
ee [el 一 人 @ 一 类 型 N 心 电 数 据 
如 下 5 ， 让 一 一 日- 类 型 J 心 电 数据 
”02 上 zg“ 
a) 对 心 电 数 据 进行 ESMD 分 解 ， 得 到 一 系列 IMF 和 一 个 趋 
全 :直下 
势 项 7 。 
ol 
b) 计算 IMF 分 量 和 7 的 能 量 已 E, 9 P1 p> ps pa ps pe pz pr 


能 量 向 量 


5 类 型 N 心 电 数据 与 类 型 ] 心 电 数 据 ESMD 能 量 向 量 图 对 比 
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2.2 ” 非 线 性 特征 

于 线性 特征 无 法 获取 信号 中 的 细微 变化 和 不 同 参数 间 的 
非 线 性 信息 ， 所 以 本 文 提取 心 电 数据 的 HOS 特征 来 挖掘 心 电 数 
据 中 的 丁 波 信息 ”…。 本 文通 过 提取 心 电 数据 的 三 阶 矩 ， 即 三 阶 
累积 量 作 为 特征 ， 计 算 公 式 如 下 : 


C7 =m(i,j) = Elx(n)x(n+i)x(n+j)] (5) 


其 中 : E[:] 表 示 统 计 期 望 运 算 符 ，1 和 j 为 时 延 参数 。 


图 6 所 示 为 正常 心 电 数 据 的 三 阶 累 积 量 及 其 等 高 线 ， 图 7 
所 示 为 含 J 波 心 电 数据 的 三 阶 累积 量 及 其 等 高 线 。 从 图 6 和 7 
中 可 对 比 HOS 可 对 两 种 心 电 数 据 进行 明显 区 分 。 本 文通 过 提取 
两 种 心 电 数据 的 三 阶 累积 量 特征 ， 得 到 40 个 特征 向 量 作为 非 
线性 特征 。 


Cum3(m,n) 


~ 20 。 
0 0 0 0 和 


延迟 量 n 


(a) 三 阶 累积 量 (b) 等 高 线 


轩 6 ”类 型 N 心 电 数据 的 三 阶 累积 量 及 其 等 高 线 


延迟 量 n 


惠及 到 mm 延迟 量 m 


(a) 三 阶 累积 量 (b) 等 高 线 


图 7 含 J 波 心 电 信号 的 三 阶 累积 量 及 其 等 高 线 


2.3 ”PCA 降 维 
为 降低 计算 量 ， 本 文采 用 PCA 来 进行 特征 降 维 。 


其 基本 思 


想 是 将 高 维 空间 的 初始 特征 经 过 变换 矩阵 ， 使 其 降 到 低 维 空间 
成 为 新 的 特征 向 量 ， 减 少 特征 维度 ， 但 仍 能 保留 初始 特征 的 主 
要 信 息 [12] 。 


将 心 电 数据 通过 ESMD 分 解 后 得 到 的 IMF 与 7 的 能 量 特征 
k 8 个 线性 特征 ， 与 心 电 数 据 的 三 阶 累积 量 共计 40 个 非 线性 
村 征 进行 融合 , 得 到 48 个 特征 向 量 , 采用 PCA 降 维 后 , 取 前 12 
个 主 成 分 作为 特征 训练 测试 分 类 器 。 


月 


ChinaXiv 合 
i 上 线性 NASM 


期 刊 > 


毋 凡 铭 ， YY SS 动 识 


2.4 ABC-SVM 
SVM 对 于 小 样本 条 件 由 很 好 的 泛 化 能 


， 适 用 于 本 文 J 波 


心 电 数据 的 识别 中 。SVM 的 目标 是 找到 一 个 最 优 超 平面 ， 能 够 
对 类 型 N 和 类 型 ] 的 心 电 数据 有 效 分 类 ， 其 表达 式 为 
y= P(X)+b (6) 


其 中 : 为 权 向 量 
进行 求解 ， 引 进 


; 卫 为 冰 值 。 找 最 优 超 平面 就 是 要 对 必 和 忆 
标 函 数 并 加 入 惩罚 项 后 为 


ol + 
i=l 
st. yO Dx)+b)>1-6 
E>0,i=12,...,n 


[Ey 


min J(@,¢)= 


| 


其 中 :多 为 松弛 


因子 ， C 为 惩罚 因子 。 


引进 Lagrange 函数 ， 


转换 为 凸 二 次 规划 寻 优 的 对 偶 问题 。 


m m 


mda- Ye, yj;(Xi, Xj), 


A (8) 
st. 二 <C 
其 中 Q > 0 为 Lagrange 乘 子 。 
对 式 〈8) 进行 求解 得 到 Ci ， 那 么 中 的 计算 公式 为 
0O= >》 wyp(0o) D(X) (9) 


最 后 ，SVM 分 类 函数 为 


Jo0 =sen(oyPx) DW)+h) (0) 


由 于 径 向 基 核 函数 仅 有 核 函 数 宽度 个 参数 需要 进行 
优化 ， 所 以 选择 其 作为 SVM 的 核 类 函数 。 a 


如 下 : 


K(x:x)=exp(—|x-—x| /oo) GD 


本 文采 用 ABC 算法 优化 与 分 类 器 性 能 相关 的 参数 C 和 ， 
建立 基于 SVM 的 丁 波 自动 识别 模型 。 具 体 步 又 为 ”: 

a) 确定 人 工蜂 群 的 群体 数量 ， 同 时 确定 侦查 蜂 所 占 比例 。 

b) 在 参数 取 值 范 围 内 随机 选择 A 个 (C ，7Y) 的 组 合 ， 
也 就 是 确定 该 A 个 组 合 即 为 蜜源 , 并 使 所 有 蜜蜂 均 成 为 侦查 蜂 。 

c) 计算 4 个 蜜源 的 品质 。 若 品质 低 , 则 侦查 蜂 不 作 变 化 继 
续 搜索 或 作为 跟随 蜂 去 高 品质 蜜源 采 蜜 ， 跟随 蜂 则 依据 蜜源 品 


录用 稿 
质 按照 一 定 比 例 去 蜜源 进行 采 密 。 
d) 当 跟 随 蜂 抵达 蜜源 ， 通 过 对 其 附近 连续 区 域 进行 采 蜜 ， 
从 而 确定 最 优点 。 
e) 所 有 的 蜜蜂 返回 蜂巢 , 计算 各 蜜源 的 品质 来 重新 分 配 各 
人 重复 上 述 a)~e), 最 终 确定 最 优 蜜源 点 , 即 最 优 的 ( C ， 
J 
3 ”结果 与 分 析 
3.1 评价 指标 
本 文通 过 灵敏 度 (sensitivity，Se) 、 特 异性 (specific, Sp) 
和 准确 率 (Accuracy，Ac) 来 评价 丁 波 自动 识别 方法 的 性 能 。 


型 N 的 比率 。 
二 3.2 结果 分 析 
之 本 文 的 数据 有 
共 约 10 400 个 ， 
动 识别 方法 进行 分 析 . 将 数据 
取 其 中 的 4 个 子 集 作为 训练 集 


个 作为 测试 集 (2 
的 一 个 样本 作为 测试 样本 ， 计 算 其 
复 5 次 , 取 5 次 的 3 
ESMD+ABC-SVM、HOS+ABC-SVM、 融 合 特征 +SVM 和 本 文 
提出 的 方法 , 各 识别 方法 中 SVM 参数 (C ,7 ) 如 表 1 所 示 ， 
平均 性 能 结果 如 表 2 所 示 。 


080 个 心 拍 ) 测 


N 


如 


均值 来 计 


估 分 类 器 的 性 


Se 表示 准确 识别 心 电 数据 类 型 J 的 比率 ，Sp 表示 准确 识别 心 


有 数据 类 型 N 的 比率 ，Ac 表示 准确 识别 心 电 数 据 类 型 J 和 类 


攻 包 含 80 个 心 电 数 据 ， 每 个 记录 持续 2 min， 
心 拍 。 通 过 五 倍 交叉 验证 法 对 本 文 提出 的 了 波 自 
平均 分 为 5 个 互 不 相交 的 子 外 
(8 320 个 心 拍 ) 训练 分 类 器 ，1 
试 分 类 器 。 每 一 次 实验 取 其 中 
性 能 评价 指标 。 这 一 过 程 重 
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能 。 分 别 对 比 


平 


ESMD+ABC-SVM 的 Se、Sp 和 Ac 平均 


和 
Se 
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Clark 等 人生 汪 厅 型 


人 台 b 
HE 
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效 
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识别 方法 ESMD+ABC-SVM 的 5 倍 交叉 验证 Se、Sp 和 Ac 
均值 分 别 为 92.7%、90.1% 和 93.8%，HOS+ABC-SVM 比 
值 分 别提 高 2.9%、5.3% 
2.5%， 将 ESMD 与 HOS 融合 ， 采 用 ABC-SVM 进行 识别 ， 
、Sp 和 Ac 平均 值 达到 最 高 , 分 别 为 96.7%、96.8%、97.3%。 
过 对 比 融 合 特 征 +SVM， 采 用 ABC-SVM 比 SVM 的 Se、Sp 
Ac 平均 值 分 别提 高 2.2%、3.2% 和 2.1%。 
步 对 比 了 目前 已 有 的 其 他 两 种 本 波 检测 方法 。 
RW 通过 形态 学 和 几何 学 特征 来 检测 本 波 ， 进 
够 对 两 种 类 型 的 J 波 进行 识别 检测 ,存在 一 定 的 局 限 性 .Wang 
人 车床 找到 引用 源 ' 提 出 一 种 利用 全 CA 和 Jscore 从 12 导 联 心电图 
自动 检测 J 波 的 方法 ， 但 算法 复杂 度 较 高 ， 且 测试 组 数据 的 
仿 结 果 不 理想 。 

相 比 于 以 上 两 种 方法 ， 本 文 从 线性 与 非 线性 角度 对 心 电 数 
提取 特征 ， 特 征 信息 较为 全 面 。 比 较 不 同方 法 对 丁 波 的 识别 
果 ， 结 果 如 表 3 所 示 。 由 表 中 数据 可 知 ， 本 文 方法 要 优 于 另 
两 种 方法 ， 对 村 波 的 识别 效果 最 好 。 

表 3 三 种 本 波 识 别 算法 的 性 能 比较 /% 


本 文 还 进 


T 


T 


文献 方法 Se Sp 


Clark et al. [15] 


Wang et al. [16] 


几何 特 


征 


NN 
区 
出 


90.5 96.5 


fPCA+J score 
融合 特征 +ABC- 


89 86 


本 文 96.7 96.8 


SVM 


表 1 各 种 识别 了 波 方法 的 SVM 参数 
方法 C 4 
ESMD+ABC-SVM 268.09 23.25 
HOS+ABC-SVM 172.43 17.02 
融合 特征 +SVM 746.26 82.18 
本 文 方法 98.85 9.73 
表 2 各 种 识别 J 波 方法 的 平均 性 能 
方法 (%) 1 2 3 4 5 均值 
Se 93.6 923 933 924 918 92.7 
ESMD+ABC- 
Sp 90.1 89.6 909 894 90.7 90.1 
SVM 
Ac 944 95.2 93.1 93.8 925 93.8 
Se 949 942 96.2 962 96.6 95.6 
HOS+ABC- 
Sp 953 95.3 95.9 96.6 94.1 95.4 
SVM 
Ac 97.1 95.9 95.3 97.7 95.4 96.3 
Se 950 953 93.3 946 94.1 94.5 
融合 特征 
Sp 94.5 94.5 93.6 925 929 93.6 
+SVM 
Ac 95.5 94.2 96.5 94.1 95.5 95.2 
Se 95.8 96.7 963 960 98.9 96.7 
本 文 方法 Sp 97.5 96.7 97.1 97.9 95.0 96.8 
Ac 98.2 97.7 96.5 966 97.6 97.3 


征 


了 一 种 丁 波 自动 识别 技术 ， 全 面 提取 心 电 数 据 特 
， 通 过 融合 极点 对 称 模 态 分 解 后 的 能 量 特征 与 非 线性 的 高 阶 


未 


中 


C 


计量 特征 并 进行 PCA 降 维 , 输入 到 ABC 算法 优化 后 的 SVM 
进行 识别 。 实 验 结 果 表 明 ， 该 方法 可 以 准确 高 效 地 识别 出 心 
数据 中 的 J 波 ， 为 医生 诊断 心脏 疾病 提供 辅助 。 
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